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ABSTRACT 

This  paper  develops  optimal  harvesting  policies  for  a  population  whose 
evolution  is  governed  by  a  single  autonomous  nonlinear  differential  equation. 
The  objective  functional  is  not  assumed  to  be  convex.  The  results  are  used  to 
discuss  optimal  policies  for  age-structured  populations  harvested  with  effort 
independent  of  age. 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  develops  optimal  policies  for  harvesting  biological 
populations.  The  age  structure  of  the  population  is  taken  into  account.  The 
policies  generally  involve  harvesting  at  maximal  effort,  not  harvesting,  and 
harvesting  to  keep  the  growth  rate  of  the  population  a  maximum. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 


summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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ON  OPTIMAL  HARVESTING  WITH  AN  APPLICATION  TO 
AGE -STRUCTURED  POPULATIONS 

★  *  + 

Morton  E.  Gurtin  and  Lea  F.  Murphy 

1 .  Introduction. 

In  this  paper  we  develop  optimal-harvesting  policies  for  a  population  whose  evolution 
is  governed  by  a  single  autonomous  nonlinear  differential  equation.1  We  present  a 

o 

formulation  of  the  infinite-horizon  problem  which  does  not  involve  a  discount  rate  and# 
using  straightforward  arguments,  are  able  to  establish  the  following  results: 

(a)  If  the  effort  is  unconstrained,  the  optimal  policy  is  to  reach  a  certain  value 
P*  of  total  population  as  quickly  as  possible,  and  then  to  hold  the  population  at  that 
value  for  all  subsequent  time.  As  one  would  expect,  P*  is  the  population  size  that 
maximizes  the  growth  rate  and  hence  corresponds  to  the  maximum  sustainable  yield. 

(b)  This  policy  is  also  optimal  when  the  effort  is  constrained  to  be  less  than  a 
constant  E,  as  lonq  as  E  is  sufficiently  large.  For  E  small,  however,  it  may  not  be 
possible  to  hold  the  population  P(t)  constant  at  P*.  In  this  case,  the  optimal  policy 
involves  a  "bang-bang"  effort  function  and  the  total  population  approaches  a  limit  larger 
than  P* . 

We  prove  these  results  directly,  without  the  aid  of  classical  control  theory.  Indeed, 
the  classical  theory  is  not  well  suited  to  this  problem,  since  we  deal,  in  part,  with 
unbounded  controls,  and  since  our  formulation  of  the  infinite  horizon  problem  is 
nonstandard. 
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This  problem  has  been  given  considerable  attention  in  the  literature  (cf.,  e.g. ,  Clark  and 
Munro  Ml,  Clark  [21,  Spence  and  Starrett  [31).  Our  results  (b)  for  and  E  small 

appear  to  be  new. 

2 

We  use  the  overtaking  criterion  of  optimality. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-004 1  and  by  the  National 
Science  Foundation  under  Grant  No.  MCS78-0 193S. 


We  use  our  results  to  discuss  optimal  policies  for  age-structured  populations 
harvested  with  effort  independent  of  age.  Using  a  simplified  form,  due  to  Coleman  and 


Simmes  [41,  of  a  general  nonlinear  model  introduced  by  Gurtin  and  MacCamy  [5],  we  are  able 
to  show  that  for  a  population  whose  age  distribution  is  initially  persistent1  -  an 
assumption  justified  for  a  population  which  has  been  evolving  over  a  long  period  of  time  - 
the  optimal-harvesting  problem  reduces  to  the  problem  discussed  above. 


2.  Age- independent  theory, 
a.  The  optimal-harvesting  problem. 

We  consider  a  species  whose  total  population 
by  a  differential  equation  of  the  form 


Figure  1.  The  growth  rate  y. 


P(t),  when  not  harvested,  is  governed 

P  =  y(P) 

We  assume  that  the  growth  rate  y 
is  a  C1  function  on  [0,«>) 
with  a 

strict  qlobal  maximum 

(1) 

at  >  0. 

We  assume  further  that,  on  (0,P^), 
y  is  strictly-positive^  and  has 
a  finite  number  of  critical 
points . 


1 

We  use  the  term  "persistent  age  distribution"  for  what  in  the  literature  is  usually 
referred  to  as  a  "stable  age  distribution."  In  a  future  paper  [61  we  will  investigate  the 
cons orm o nces  of  arbitrary  initial  conditions. 

This  assumption  can  be  weakened  to  include  "critical  depensation." 

»<5sumo  that,  when  P ^  <  P„,  y  >  n  on  {Pry?^* 


Indeed,  we  need  only 


Suppose  that  at  t  »  0  the  population  has  in.'.tial-value 

P(0>  =  PQ  >  0  ,  (2) 

and  that  subsequently  individuals  are  harvested  at  a  rate 

E(t)P(t) 

with  E  a  nonnegative  function  representing  the  (harvesting)  effort.  Then 

P  =  y(P)  -  EP  ,  (3) 

and  using  this  equation  and  (2)  the  yield 

/q  E(t)P(t)dt 

on  a  time  interval  (0,T)  can  be  written  as  a  functional  of  P: 

Yt!P)  -  PQ  -  P(T')  +  J*  y(P(t) )dt  .  (4) 

Let  us  agree  to  use  the  term  path  for  a  strictly-positive,  piecewise  continuous, 
right-continuous,1  oiecewise  C1  function  P  on  [0,«*).  Assume  that  we  are  given  a 

family  A  of  paths.  Then  an  optimal  path  relative  to  A  is  s  path  P  e  A  with  the 

.  * 

follow'-q  property:  given  any  P  e  A,  P  P  , 

Yt(P*>  >  Yt(P)  (5) 

for  all  sufficiently  large  T  (that  is,  for  all  T  larger  than  some  T^  *>  T^(P)).  Note 
that,  since  (S)  is  strict,  there  is  at  most  one  optimal  path. 

2 

There  are  many  Dossible  objective  (profit)  functions  with  which  one  can  work.  For 
convenience,  we  have  chose  to  optimize  the  total  vield.  Our  analysis  goes  through,  almost 
without  change,  when  a  cost  term  proportional  to 

1  *’”* 

T.e.,  p(t)  =  P(t  ).  This  assumption  removes  ambiguity  at  points  of  discontinuity. 

o 

A  thorugh  discussion  of  objective  functions  aooropi iate  to  harvesting  problems  is  qiven  by 
Clark  r 2 1 . 
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1 


fo  E(t)dt 


is  subtracted  from  (4). 

Remark.  The  inf inte-horizon  problem  is  usually  formulated  with  the  aid  of  a  discount 

-  ift 

factor  e  (k  >  0);  one  then  maximizes  the  total  yield 

f“  E(t)P(t)e_<tdt  . 

This  formulation  renders  the  yield  in  the  immediate  future  more  important  than  that  of  the 
distant  future.  To  the  contrary,  our  procedure  models  situations  in  which  one  is  willing 
to  accept  a  sub-optimal  yield  initially,  in  order  to  eventually  produce  the  best  yield 
possible.  The  use  of  a  discount  factor  may  be  sound  from  an  economic  point  of  view,  but 
our  policy  is  certainly  more  appropriate  to  a  long  range  conservation  program.  It  is 
interesting  to  note  that  discounting  generally  leads  to  a  smaller  size  for  the  ultimate 
population. 

b.  Optimal  harvesting  when  the  effort  la  allowed  to  be  unbounded. 

Our  first  step  will  be  to  define  the  class  A  of  paths.  To  begin  with,  we  require 
that  the  effort  be  nonnegative;  in  view  of  (3),  this  is  insured  by  the  constraint; 

( A  ^ )  for  any  t  >  0  at  which  P  is  differentiable, 

P(t)  <  y<P(t)  )  • 

Next,  note  that  YT(P)  is  well  defined  for  all  paths  P,  even  those  that  suffer  jump 
discontinuities.  If  individuals  are  harvested,  rather  than  stocked,  the  population  after 
such  a  jump  must  be  less  than  the  population  before.  We  therefore  add  the  second 
constra lnt : 

(A^)  P(0)  <  Pg,  and  for  any  t  >  0  at  which  P  jumps, 

P(t)  <  P(t")  . 

To  interpret  the  differential  equation  (3)  on  paths  with  jump  discontinuities,  we 


rewrite  if  in  the  form 


(  in  P) 


^  y(p>  -  e  . 


For  this  equation  to  make  sense  for  functions  with  jumps,  E  must  be  interpreted  as  a 
distribution.  Indeed,  if  t^  (k  =  1 , 2, . . . , n)  denote  the  times  of  discontinuity  of  P, 
then 


E(t> 


E  <t)  +  <5(  t )  in 


P(0) 


+  l  5(  t-t  )  in 
k  k 


P(t.  ) 


with  Eg  piecewise  continuous  and  6  the  delta  distribution.  Note  that,  with  E  defined 
in  this  manner,  (Aj)  follows  from  the  requirement  E  >  0. 

Since  y  has  a  strict  maximum  at  PM,  we  would  expect  the  optimal  policy  to  involve 
reaching  PM  as  quickly  as  possible  and  then  remaining  there  for  all  subsequent  time.  For 
Pg  >  PM  this  is  easily  accomplished  by  jumping  immediately  to  PM,  since  such  a  jump  is 
consistent  with  the  constraint  (Aj).  For  PQ  <  PM  a  jump  from  PQ  to  PM  is  not 
possible)  as  we  shall  see,  here  the  optimal  way  to  reach  PM  is  to  refrain  from 
harvesting.  With  this  in  mind  we  introduce  the  following  definition:  the  f ree-qrowth 
curve  starting  at  ZQ  e  ( 0 , PM )  at  time  tQ  >  0  is  the  solution  Z(t),  t  >  tQ,  of  the 
initial-value  problem 

z  =*  y(Z)  , 

Z(tg)  =  ZQ  . 

Note  that,  since  y  >  0  on  (0#PMJ,  Z  is  strictly-increasinq  and  reaches  PM  in  finite 
time. 

Theorem. 1  Relative  to 

A=  {paths  consistent  with  ( A  ^ )  and  ( A^ )  } 
the  optimal  path  is  given  by 


1 

Cf.,  e.q.,  Clark  and  Munro  Ml,  Clark  [2],  Chapter  2. 
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t  >  0 


(6) 


p*(t>  *  PM, 

for  >  P  -  and  by 

-  0  M  - *- 


P  (t)  = 


Z(t),  o  S  t  <  t 


p  ,  t  <  t  < 
.  M  M 


(7) 


for  P(1  <  P^.  Here  Z  is  the  free^growth  curve  starting  at  Pq  at  time  t  -  0,  while 
is  the  time  at  which  Z  reaches  Pu. 

“  '■  '  —  ■  ■  —  ■  ■  -  —  ■  M 


I'l.im'i:  timal  policies 


By  (3),  the  effort  E  required 
to  hold  P(t)  constant  at  PM  is 
given  by 

e  -  y(PM)/PM. 

Thus  for  the  path  (6), 

P 

E(t)  =  y(P„)/P„  +  6(tUn  —  , 

M  M  PM 

while  for  (71, 


E(t> 


/*•  0‘l<t" 

l  v<V'V 


f'ur  v-roof  of  the  theorem  begins  with 
Assertion  1.  Let  P  6*  A  be  optimal.  Then: 

f  i  )  p(  M  k  r  for  all  t  >  0 ; 

(  i  i  >  p(t)  Ps<  for  all  sufficiently  larqe  t. 

Proof,  Assume  that  (  i)  is  not  true,  so  that  P  >  P^  on  some  set  J  c  [0,«>).  Define 

o  p  Ah*/  -  min[P  , P(t)}.  Then,  since  G  <  P,  it  follows  from  (1)  that 

M 

Y  (T)  -  Y  (P)  =  '  [y(P  )  -  y(P( t) )1  dt  +  P(T~)  -  G(T_)  >  0 

T  T  fO,T)  J  ” 


whenever  T  >  inf  J,  which  contradicts  our  hypothesis  that  P  be  optimal.  This  proves  (i). 
Before  beginning  the  proof  of  (ii)  note  that,  by  (1)  there  exists  an  L  e  (0,P„)  such 

that 

p,  <  p2,  p2  e  ->  y(P,>  <  y(P2)  •  (8) 

Our  first  step  in  establishing  (ii)  will  be  to  show  that 

lim  P(t)  -  Pu  .  (9) 


Assume,  to  the  contrary,  that  (9)  is  not  true.  Then,  by  (i), 

lim  inf  P(t)  <  , 

M 


and  there  exists  a  K  e  (L,P^)  and  an  increasing  sequence 


{t  }  with  t 
n  n 


and 


P(tn)  <  K 

for  all  n.  Let  Z  denote  the  free-growth  curve  starting  at  K  at  t  =•  0,  let  x  denote 
the  time  at  which  Z  reaches  PM,  and  let 


Zn(t)  *  «*-V  ' 

so  that  z  is  the  free~growth  curve  starting  at  K  at  time  t„ . 


Further  assume,  without  loss  in  generality 

that  tn+i  >  tn  +  x  for  all  n.  (If 

{t  }  does  not  have  this  property,  some 
n 

subsequence  does.)  Since 

P(t  )  <  Z  (t  )  , 

n  n  n 

Z  =  y ( Z  ),  P  <  y(P)  , 
n  n 


Figure  3.  Bounding  p  by  Z  . 

n 


and  since  the  jumps  in  P  are  downward, 
we  may  conclude  from  a  well  known  comparison 
theorem^  for  differential  inequalities  that 


1 

Cf.,  e . q. ,  Hartman  [7],  Thm.  4.1,  p.  26. 
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(10) 


P  <  Z  on  ( t  ,  t  +  <1 
n  n  n 

* 

Thus,  for  P  defined  by  (6)  or  (7)  and  t  sufficiently  large. 


j\n+1  [y(P* (t) )  -  y(P(t))]dt  =  f  n+1  !y(  P  )  -  y(P(t)))dt 
t  t  m 


>  r  tv(p„)  -  V(z  (t)))dt  =  f*  ty (P  )  -y( Z( t) ) ) dt  =  D  >  0 

t  “  n  ■  i)  *1 


where  we  have  used  (8);  hence 


F(T)  =  Jg  [y(P  (t))  -  y(P(t)>ldt  -►  +»  as  T  + 


Thus,  since  P  >  0, 

*  _ 

Y  (P  )  -  Y-(P)  =  F(T)  +  P(T  )  -  p  ♦  +oo  as  T  ♦  «,  , 

T  T  M 

which  again  contradicts  the  assumed  optimality  of  P.  Thus  P  satisfies  (9). 
Next,  by  (8)  and  (i),  there  exists  a  time  a  for  which 


Let  Z  be  the  f r ee-gr^wt h  curve 


starting  at  P(a)  at  tine  t  =  a, 


let  <  denote  ♦’he  :\me  at  which  Z 


reaches  PM,  and  let  G  p  4  be 


t  defined  by 


Figure  4.  Construction  of  G. 


P(  t )  ,  0  v,  t  <  a 
G(  t)  =  (  Z(  t ) ,  a  <  t  <  * 


V  *  *  1  < 


Yt(G)  -  Yt(P)  =  fa  !y(0(t) )  -  y(P(t))]dt  +  P(T_i  - 


por  all  T  >  *r.  Further,  the  argument  leading  to  (10)  here  tells  us  that 


P  <  Z  on  [a ,  <] 


Suppose  that  (ii)  is  not  satisfied-  Then,  by  (8),  (11),  and  the  properties  of  G, 
y(G( t ) )  >  y ( P ( t ) )  for  all  t  e  la,*)  , 

the  inequality  beinq  strict  on  an  open  set.  Thus,  since  P(t)  -*■  P  as  t  «,  {  '  ■'  mu 
be  strictly  positive  for  T  sufficiently  large.  Since  this  cannot  be  so,  (ii)  must  be 
valid. 

In  view  of  Assertion  1,  it  suffices  to  establish  optimality  within  the  class 

A  =  {P  e  A  |  P  <  PM  on  [0 ,«) ,  P( t )  *  Pu 
MM  M 

for  all  sufficiently  large  t) 

*  * 

Assume  that  >  P  .  Let  P  be  given  by  16)  a nd  choose  P  8  A*,,  P  ¥  P  .  Then 

P(T)  *  P*(T)  *  PM  for  all  sufficiently  2arge  T,  and  for  all  such  T, 

yt<p*1  *  VP)  =  .>1  [y(pm’  "  y(P(t) ))dr  >  0  , 

* 

hence  P  is  optimal. 

To  complete  the  proof  we  must  establish  the  optimality  of  (7)  when 

P0  <  PM  • 

With  this  in  mind,  we  introduce  the  following  notation.  Let  [A,C]  (0,P^},  let  Z 

denote  the  free-growth  curve  starting  at  A  at  t  -  0,  let  *'  denote  the  time  at  w''ic 
7  reaches  C ,  and  let  [a,*?1  [h,°°)  with 


- 


c  >  a  +  k 


M2) 


Further,  let  H  be  defined  on  ta,c)  by 


H(t) 


( 


Z(t-a),  a  <  t  <  a  +  < 


a  +  >c  <;  t  <  c 


if  Y(C)  >  y(A); 


H  ( t )  = 


A,  a  <  t  <  c  -  < 

Z(t-c+<),  c  -  <  <  t  <  c 


(13) 


li.iuro  \> .  optimal  transitions.  y(C)  <  Y(A).  We  call  H  the 

optimal  transition1  from  A  to  c 
during  [a , cl ;  <  is 

the  transition  time  of  H,  c  -  a  -  <  the  rest  time  of  H.  Finally,  a  path  P  e  AM  is  a 

r  i  so  during  [a,cf  if 

P( a)  <  P(t)  <  P(c)  for  all  t  e  (a,c)  .  (14) 

Assert  ion  2.  Let  P  f*  Am  be  a  rise  during  (a,c)  «■  Then: 

(i)  There  is  an  optimal  transition  H  between  P(a)  and  Pfc)  during  [afc]. 
(ii)  If  the  rest  tine  of  H  vanishes,  then  P  *  H  on  (a, cl. 

!iii)  If  y  in  monotone  on  lP(a),P(c)1,  then 

y ( H ( t ) )  >  y(P(t))  for  all  t  e  [a, cl  , 
vo.  t.h  strict  inequality  on  an  open  set  if  the  rest  time  of  H  is  nonzero. 
Proof,  w  e  will  rive  the  proof  on)  v  for 

y ( P( c) )  >  v(P( a) ) 

Pi ea r  1  y , 

P(t.)  <  7  ( t  -  a ) ,  a«;t<a  +  K  (is) 


Pf  ,  ,  Spence  anil  Jltarrett  f3l,  n.  3^3, 
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(cf.  (10)).  Thu9  c  is  consistent  with  (12)  and  the  optimal  transition  (13)  is  well 
defined.  Assume  (for  the  remainder  of  this  paragraph)  that  the  transition  time  <  of  H 
equals  c  -  a.  Then  H(t)  *  z(t-a)  for  all  t  €  (a,c),  and  since 

H(c)  -  P(c)  , 

H  -  y(H),  P  <  y(P)  , 

it  follows1  that  P  >  H  on  (a,c);  this  inequality  and  (15)  yield  P  *  H  on  (a,c]. 

To  prove  (iii)  let  y  be  monotone  and  hence  strictly  increasing  on  [P(a),P(c)]. 

(The  monotonicity  is  strict  since  y  has  at  most  a  finite  number  of  critical  points  in 
(0,P„).)  It  is  clear  from  (15)  and  (14)  that  P  <  H  on  [a,c) ,  and  if  the  rest  time  of 
H  is  nonzero,  that  P(a+<)  <  H(a+<)  =  P(c).  This  yields  the  validity  of  (iii). 

We  are  now  in  a  position  to  establish  the  optimality  within  AM  and  hence  A  -  of 
the  path  P*  defined  by  (7).  Thus  choose  P  e  AM,  P  P*,  and  let  TM  be  the  least  time 
for  which 

P(t)  =  Pw  for  all  t  >  T  .  (16) 

M  M 

Since  all  jumps  in  P  are  downward,  P  must  take  on  all  values  in  tPQ,PMl.  Let 
C1'C2 . Cn  with  P0  “  C1  <  C2  <•••<  Cn  =  PM 


I 


Flour*?  6.  Decomposition  of  p. 


r*f.,  e.ri.,  Hartman  (71,  Remark  1,  p.  76. 


1 


denote  the  extreme  points  of 


n)  denote  the  last 


y  in  [P0,PM1.  let  <k  =  1,2 
time  in  [0,TM1  at  which  P  takes  on  the  value  Ck  (note  that  £ =  T^),  let  tk 
(k  =  2, 3, .... n)  denote  the  first  time  after  1  at  which  P  has  the  value  Ck,  and 
define  r  =  0.  Then 

P  is  a  rise  during  ( i,  ,  T,  .]  (k  =  1,2,...,n-1)  .  (17) 

k  k+1 

Consider  the  path  R  e  AM  obtained  from  P  by  replacing  each  rise  (17)  by  the 

corresponding  optimal  transition.  Then,  since  y  is  monotone  on  each  [P(£  ),  P(lk+^)1» 

we  may  conclude  from  Assertion  2  that 

y( R( t ) )  >  y(P(t))  for  all  t  >  0 

(Note  that  P  and  R  coincide  on  each  !  T,  ,  l.  ]  • )  Thus 

k  k 

Y  (R)  >  Y  (P)  for  all  T  >  T  .  (18) 

T  T  M 

We  shall  complete  the  proof  by  showing  that 

Yt(P*)  >  YT(P>  for  all  T  >  T^  .  (18) 

*  » 

For  convenience,  we  consider  separately  the  two  cases:  R  =  P  ,  R  r  P  . 

Case  1  (R  »  P*).  Here 


Tk  =  lk  ~  1,2, ...,n)  and  the  rest  time  - 
for  each  rise  except  the  last  -  vanishes 


(20) 


Thus,  by  (ii)  of  Assertion  2,  R  and  P  coincide  outside  of  the  last  rise-interval 

[£  ,r  ).  The  rest  time  associated  with  this  last  rise  cannot  vanish,  because  if  it  did, 

n- 1  n 

* 

then  R  and  P  would  coincide  everywhere,  an  impossibility,  since  R  =  P  *  P. 

Consequently,  we  may  conclude  from  (iii)  of  Assertion  2  and  (8)  that 

Y(P)  <  Y  ( R )  (=  Y  (P  ))  for  all  T  >  T  . 

T  T  T  M 

Case  2  ( F  ^  P*).  Let  T  >  Tu  =  £  .  By  (7),  the  integral 

-  “  n 

(p  y(P  (t) )rtt 
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* 

is  equal  to  (T-T^)y(P^)  plu'  the  inteqral  of  y(P  )  over  the  free-qrowth  segment  of 

p*.  Since  the  free-growth  equation  is  autonomous,  this  latter  integral  equals  the  integral 

of  y(R)  over  the  free-growth  portions  of  R.  Thus 

Yt(p'')  -  Yt(R)  =  (T-TM)y(P„)  -  It(R>  ,  (21) 

where  IT(R)  is  the  integral  of  y(R)  over  the  set  SI  consisting  of  the  rest  portions  of 

the  optimal  transitions  that  comprise  R  and  the  union  of  the  intervals  (t  ,£^), 

(t  ,£  ,T).  On  SI,  y(R)  <  y(P  )  and  (as  we  shall  show  at  the  end  of  the  proof) 

2  2  n  M 

SI  contains  an  open  set  0  with 

y(R(tl)  <  y(P„)  for  all  t  e  0  .  (22) 

Thus,  since  the  measure  of  ft  is  T  -  TM, 

It(R)  <  (T-TM)y(PM)  , 

and  (21)  yields  YT(P*)  >  YT(R).  In  view  of  (IB),  this  inequality  implies  (19). 

To  verify  (22)  note  first  that  (20)  cannot  be  satisfied,  for  if  it  were,  then  R 

would  equal  P* .  Thus  at  least  one  interval  ( , £^1  or  at  least  one  rest  interval  with 

R  /  p„  has  nonzero  duration.  This  interval  clearly  contains  an  open  set  C  on  which  (22) 

is  satisfied.  (If  the  interval  in  question  is  tr  ,1  ],  the  existence  of  0  is  insured 

n  n 

by  the  equality  of  R  and  P  on  an(^  the  ^act  that  is  the  smallest 

time  consistent  with  (16)). 

Remark  1.  The  construction  P  -  obtained  from  P  by  replacing  each  rise  (between 
extrema  of  y)  by  its  corresponding  optimal  transition  -  can  be  defined  (in  an  obvious 
manner)  for  anv  function  P  with  domain  an  interval  fa,h),  as  Inna  as  P(a)  <  P(b)  and 
P(t)  <  P(b)  for  all  t  f?  (a,b).  We  will  refer  to  this  const  ruct  ion  as  the  qua  si -opt  ini?,  or 
of  p  on  fa ,  b’ . 
c .  Bounded  effort. 

The  problem  becomes  more  complex  when  the  effort  is  restricted  by  ar,  upper  bound  F. 

We  require  that  E(t)  <.  F  for  all  t,  or  equivalently, 

( A  ^ )  P(t )  >  y(P( t ) )  -  EP(t) 

A  path  which  obeys  (A^)  and  ( )  is  necessaril v  continuous.  We  therefore  lismiss  (AO  an  1 
restrict  our  attention  to  paths  consistent  with  (Aj),  (A-j),  and  the  initial  condition  (2). 


In  addition  to  assuming  that  y  is  a  function  on  [0,«)  with  a  finite  number  of 

critical  points,  we  assume 

lim  sup  y(P>  <  0  (23) 

P  >  oo 


Then  there  is  a  largest  value  P^  ( <  «)  for  which 

V(PS)  =  iPs  - 

Further,  by  (23),  y  restricted  to  has  a  maximum;  we  assume  that  this  maximum 

occurs  at  exactly  one  point  P  e 

M  S 

y(P  )  >  y(p)  for  all  p  >  p  ,  p  ^  P  .  (24) 

M  S  “ 

Finally,  we  assume1  that  y  >  0  on  (0,Pg>.  We  use  the  notation 

S  =  y(Ps),  M  =  y(PM)  ;  (25) 

then  M  >  S. 

In  addition  to  free  growth  curves,  the  optimal  path  involves  maximal-effort  curves; 

that  is,  solutions  X(t),  0  <  t  <  <»,  of  the  equation 

X  =  y(X>  -  EX  .  (26) 

2 

we  aqain  aopeal  to  a  comparison  theorem  to  note  that  if  P  satisfies  (r3)  and  P(tp)  « 

X ( t q  ) ,  then 

P  <  X  on  (0,tQ),  P  >  X  on  Itg/®)  •  (27) 

For  Pn  <  PM,  we  will  use  the  following  notation: 

U1  =  inf {p  >  P0  I  y(p)  >  M}  , 

and  for  each  value  of  n(n  =  1,2,...)  for  which  the  underlying  set  is  nonempty, 

V  =  inf  {p  >  U  |  y(p)  <  M)  , 
n  n 

U  =  inf  (p  >  V  |  y(p)  >  M} 
n+ 1  n 


Then 


1 

i  ; i~r*  i-'i  may  be  weakened  as  indicated  in  Footnote  2,  p.  2. 
' ' f . ,  e.g.,  Hartman  [7},  Thm.  4.1  and  Remark  1,  p.  26. 
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(28) 


y  >  M  on  (U  ,V  )  , 

n  n 

and,  by  (24),  V  (,  P  ,  for  all  n.  Since  y  has  a  finite  number  of  extrema,  it  is  clear 
n  S 

that  there  is  a  finite  number  N  (say)  of  U  '  s  and  the  same  number  of  V  's. 

'  n  n 


We  consider  two  cases. 


Case  A,  PM  is  not  a  local 


maximum  of  y.  Then,  by  (24), 


PM  =  Ps,  y  <  S  on  (P  ,«•)  , 


P0  U1  V1 


Figure  7.  The  function  y  for  Case  A. 


Theorem  A.  Relative  to 


■  p  and  VN 


A  «  (paths  consistent  with  ( A^ ) ,  { A3 ) ,  and  (2)} 


the  optimal  path  Is  P  ,  as  described  below. 


(a)  If  Pn  <  P„,  then 


j'  0  ‘  *  <  u, 

j  Xt(t),  U1  <  t  <  v1 


ZN<t)’  Vi  4  "  <  °N 


XN(t),  UN  <  t  < 


Z j  the  free-growth  curve  through  Pq  at  t  =  0, 

Zn  the  free-growth  curve  through  t  =  vn-1' 

the  time  at  which  Z „ ( u_ )  =  U  , 
n  -  n  n  n 


the  time  at  which  X„ ( v )  =  V 
- - - n  n  n 
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(b)  _!£  pQ  >  P  ,  then  P  la  the  maximal-effort  curve  through  PQ  at  t  »  0. 


We  justify  the  vn  by  noting  that, 

for  n  <  N,  y(p)  -  Ep  >  0  on  [On#Vnl , 

so  that  X_  reaches  V„  in  a  finite 
n  n 

time.  The  optimal  paths  are  graphed  in 
Figure  9.  Before  proving  the  theorem  for 
Case  A,  we  state  the  results  for 


a  vp 

N  NS 


Figure  8.  The  function  y  for  Case  B. 

Case  B.  PM  is  a  local  maximum  for  y  (Figure  8). 

Theorem  B.  Relative  to 

A”  (paths  consistent  with  ( A^ ) ,  ( A^ ) ,  and  (2)} 
the  optimal  path  is  P  ,  as  described  below 


<»>  II  po  <  V 


Z,(t),  0  <  t  <  u. 


X,(t),  u,  <  t  <  v. 


Vl(t)’  VN  <  *  < 


fcM  <  < 


Here  Zn<  Xn ,  un  and  vp  are  interpreted  as  in  Theorem  A,  ZN+1  is  the  free-qrowth  curve 
through  VN  at_  t  =  vN,  and  t*,  is  the  time  at  which  ZN+1(tM)  »  PM. 


<b)  If  P0  >  PM  >  V 


X  ( t )  ,  0  <  t  <  t 


V  < 


where  X  19  the  maximal-effort  curve  through  PQ  at^  t  -  0,  and  tM  is  the  time  at 

which  X ( )  -  PM.  (Since  X(®)  *  P  <  P  ,  t  <  ) 

-  M  "1  -  S  MM 

(c)  If.  and  Ps  =  PM,  then  P*  is  the  maximal-effort  curve  through  Pn 


at  t  ~  0 . 
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Figure  9.  Optimal  paths. 

We  will  prove  the  theorem  only  for  Case  A(a).  Since  the  technical  aspects  of  the 
proof  are  rather  complicated,  we  offer,  as  a  preliminary,  the  following  intuitive  argument 
in  support  of  the  optimality  of  the  path  (30). 

We  first  ask,  is  there  an  optimal  equilibrium  value  for  the  population?  By  (A^),  any 
equilibrium  value  p  must  satisfy  v(p)  <  Ep.  It  is  clear  from  Figure  8  that  the  "best” 
such  value  is  p  =  Pg.  We  therefore  expect,  and  will  indeed  show,  that  the  optimal  path 
asymptotically  approaches  Pg. 

We  next  ask,  what  is  the  best  way  to  make  the  population  grow  from  its  current  size  to 
the  optimal  size  P^?  We  know  that  the  optimal  asymptotic  yield  rate  is  S.  Thus  if  the 
current  yield  rate  is  less  than  S,  we  refrain  from  harvesting  iu  hasten  the  qrowth  of  the 
population  to  a  size  which  will  produce  a  higher  yield.  On  the  other  hand,  when  the 
current  yield  rate  is  higher  than  S,  we  are  in  an  unusually  profitable  situation. 

Although  the  population  will  inevitably  grow  beyond  this  profitable  stane,  we  can,  bv 

* 

harvesting  with  maximal  effort,  prolonq  the  period  of  high  profit.  The  path  P  as 
defined  in  (30)  behaves  in  the  manner  iust  described. 


-17- 


which  will  be  the  case  if 


Remark  2.  If  y  has  a  strict  global  maximum  at  P^  >  Pg, 

E  is  sufficiently  large,  then  it  is  possible  to  harvest  so  that  P  is  held  constant  at 
PM.  For  this  case  Theorem  B  is  completely  analogous  to  our  theorem  for  unbounded  effort. 
For  Pg  <  PM  the  optimal  path  is  a  free-growth  curve  followed  by  P(t)  5  P  ;  for  pg  >  PM 
the  optimal  path  is  a  maximal-effort  curve  (which  replaces  the  jump  to  PM)  followed  by 
P(t)  2  PM. 

We  now  begin  our  proof  of  Theorem  A(a).  Thus  let 

P0  <  PS  • 

Assertion  3.  Let  P  e  A  be  optimal.  Then: 

(i)  lim  P(t)  *  P  ; 
t  -*03 

(ii)  P  is  nondecreasing. 

Proof .  Note  that  P(t)  <  P  ,  0  <  t  <  ®,  since  otherwise,  the  path  Q  e  A  defined  by 

Q(t)  =  min{P(t) ,Pg} 

would  be  better.  (For  convenience,  we  use  the  terminology:  Q  is  better  than  P  if 
Y,p(C)  >  Y^fp)  for  all  sufficiently  large  T.)  Assume  first  that  P(tg)  S  tUn,VNl  at  “ome 
time  tQ.  Let  X  be  the  maximal-effort  curve  through  P(tg)  at  tQ.  By  (26),  X(t)  •*  Ps 
as  t  +  «>;  thus,  since  X  <  P  <  Pg  on  (t^,*),  ( i)  is  satisfied.  Assume  next  that 

P(  t)  <  tJ,,  for  all  t.  Lot 

Sg  =  {p  e  [0,uN)  l  y(P)  >  s  -  e)  .  ( 3 1 > 

Then  there  exists  an  e  >  0  such  that 

y(p)  >  Ep  for  all  p  e  S .  (32) 

Moreover,  since  y  has  at  most  a  finite  number  of  critical  points,  Sg  is  the  union  of  a 
finite  number  of  compact  intervals.  By  (26),  (27),  and  (32),  if  P  enters  any  such 
interval,  it  remains  only  a  finite  time,  and  it  never  returns.  Thus  P(t)  lies  outside 
for  all  sufficiently  large  t,  and  (25)1  and  (31)  yield 

y(Ps )  -  y (P( t) )  >  e  (33) 

for  all  such  t.  let 
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Q(t)  - 


0  <  t  <  t„ 

s 

t  <  t  <  » 

s 


where  Z  is  the  free-growth  curve  starting  at  Pq  at  t  »  0,  while  tg  is  the  time  at 

which  z  reaches  Pg.  Then  a  trivial  computation,  based  on  (21),  shows  that 

lim  t*T<Q>  -  Y?  ( P ) )  -  +»  , 

T**» 


so  that  Q  is  better  than  P.  This  contradicts  the  optimality  of  r>  and  the  proof  of  (1) 
is  complete. 

•  -V 

(ii)  We  will  show  that  P(t  )  >0  for  all  t.  Assume,  to  the  contrary,  that 

•  +  — 

P(t1 )  <0  at  some  time  t ^ ,  and  let  ■  P(t^).  Then  by  (A^),  y(Pj)  <  EP^ ,  and, 

since  y(Pg)  32  EPS»  pi  <  pg*  In  addition, 

q  «  sup {p  e  [O.P^  |  y(p)  =  Ep}  , 

T  »  inf {p  e  (P^Pg]  I  y(p)  “  Ep} 

exist  (recall  y(0)  >  C),  and 


y(p)  <  Ep  for  all  p  e  (q,r)  , 

0  <  q  <  P1  <  T  <  Pg  . 


(34) 


(Cf.  Figure  10.  The  inequality  t  <  Pc  is  a  consequence  of  (28)  and  (34)^). 

Let  X  be  the  maximal-effort  curve  starting  at  at  t  *  t^.  Since  the  maximal- 

effort  curve  through  q  remains  constant  at  q,  X  >  q  on  (t^,®);  hence  (27)  yields  the 
concl us  ion 

P  >  q  on  ftj,®)  .  (35) 

Let  +•  l  be  the  first  time  after  t1  at  which  P(t^+£)  -  P^.  (Since  P1  <  Pg, 

the  existence  of  H  >  0  is  assured  by  (i).)  Then 

P  <  P1  on  ( t 1 ,tl+£)  .  (36) 

By  (34)2,  (35),  and  (36),  q  <  P(t)  <  T  on  (t^#t^+t),  so  that,  by  (34)^, 
y ( P ( t )  )  <  EP(  t )  <  Et  =  v (  t )  for  all  t  e  (t^t^  +  e)  . 

Thus  the  path  Q  8  A,  pictured  in  Figure  11  and  defined  by 
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Q(  t) 


'  Q(t),  0  <  t  <  t 

P(t+£>,  t1  <  t  <  t  -  £ 

t2  -  £  <  t  <  t2 


V  P(t)  ,  t  <  t  <  *  , 


where  t2  (>  t  +£)  is  the  first  time  after  t^  at  which  P(t£)  *  T/  better  than 

•  4- 

Hence  P(t  )  <0  is  not  possible.  This  completes  the  proof  of  Assertion  3. 


Figure  10  Figure  11 


The  followinq  definition  will  be  useful.  Let  f  be  a  continuous,  nondecreasinq 
function  on  (0,«>),  and  let  F  belonq  to  the  range  of  f.  Then  f  stretched  at 

F  by  amount  i  >  0  is  the  function 
j  ?  defined  by 


f  ✓  r  f(t>'  °  ^  <  fc0 

F  _  x - ^  h t)  =  <  f,  t0  <  t  <  t0  +  i 

t  f(t-£),  t_  +  £  £  t  <  » 

y  o 

F  F 

trt  t  +£ 

0  0 

Figure  12.  f  stretched  at  F. 
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where  tQ  is  the  first  tine  at  which  f(tg)  *  F.  A  typical  stretched  function  is  shown  in 
Figure  12.  Of  course,  by  successive  applications  of  this  definition  we  may  stretch  a 
function  at  a  finite  number  of  points  in  its  range. 

* 

We  are  now  in  a  position  to  establish  the  optimality  of  the  path  P  defined  by 

* 

(30).  In  view  of  Assertion  3,  it  suffices  to  show  that  P  is  better  than  any  path 

P  e  A,  P  i*  P  ,  consistent  with  (i)  and  (ii)  of  Assertion  3.  Let  P  be  such  a  path,  let 

u  (n  *  1 , 2 , . •  • , N )  denote  the  time  at  which  P  crosses  U  ,  and  let  v  (n  =  1,2,...,N-1) 
n  n  n 

denote  the  time  at  which  P  crosses  Vn»  (These  times  are  unique*,  indeed,  if  t  denotes 

any  one  of  these  times,  then  P(t)  >  y(P(t))  -  EP(t)  >  0.)  Let  P  be  P  stretched  at 

V. , V-, . . . , V„  .  with  A  -  (v  -u  >  the  amont  of  the  stretch  at  V  .  Here  A  is  the  time 
12  N-*  I  n  n  n  n  n 

it  takes  the  maximal-effort  curve  starting  at  Un  to  reach  Vn«  Finally,  let  u^ 

(n  *  1,2,...,N)  be  the  time  at  which  P  crosses  U  ,  so  that  v  =  u  +  A 

n  n  n  n 

-►  •* 

(n  *  1,2,...,N-1)  is  the  last  time  P  takes  on  the  value  Vn.  Then,  since  P  ^  P  , 

P  t  P*  ?  (37) 

in  addition,  a  simple  analysis  shows  that 

Um  fy  (p)  -  y  (P)l  =  o  . 

T>«5 

It  therefore  suffices  to  show  that 

*  -»• 

P  is  better  than  P  (3R) 

♦ 

Consider  the  path  R  obtained  from  P  as  follows: 

(I)  On  [u  ,v  ]  replace  P  by  the  maximal-effort  curve  starting  at  U_  at  t  -  u  . 

n  n  n  n 

(Note  that  P  is  defined  so  that  this  maximal-effort  curve  reaches  V  at  v  .) 

o  n 

(IT)  On  (v  ,u  ]  replace  P  by  its  quasi-optimizer  (cf.  Remark  1).  Here  (I)  and 
n- 1  n 

■»  * 

(II)  hold  for  n  =  1,2,...,N  with  vn  =  0  '^1  v  =  «>. )  The  paths  P  and  R  are  shown 

v  n 

in  Figure  13.  Note  that  these  oaths  generally  do  not  belong  to  A. 


Figure  13.  The  curves  P  and  R. 


Assert  ion  4. 

lim  inf  [Y  (R)  -  YT<P)3  >  0 
T 


with  strict  informality  if  R  /  P. 


Proof .  Let  T  >  u^  and  define 


T  *  U  fv  .u  ]  ,  J  =  U  [u  ,v  1,  Jm  =  J  n  [0,T] 
<  n-1  n  .  n  n  T 

n~  1  n-1 


ICA)  =  .  [y ( R( t ) )  -  y(P(t))]dt  , 

•  A 


(3^ 


so  that  (0,Tl  =  I  u  JT  and 


Y^(R)  -  Y^fP)  =  P(T)  -  R(T)  +  1(1)  +  I  (  T? ) 


Ffu  )  =  R(u  ),  (?7)  implies  P(T)  >  R(T)*  Thus  to  prove  Assertion  4  it  suffices 


to  show  that  (40)  and  (41)  below  are  valid: 

1(1)  >0  with  strict  inequality  if  R  ft  P  on  I  ;  (40) 

lim  inf  /(J^)  >  0  with  strict  inequality  if  R  j*  P  on  J  .  (41) 

T  ♦  ® 

Conclusion  (40)  is  an  immediate  consequence  of  (it)  and  (iii)  of  Assertion  2  and  the 
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fact  that  R  ^  P  on  I  only  if  at  least  one  of  the  optimal  transitions  associated  with 
R  has  a  nonzero  rest  time. 

we  now  prove  (41).  On  an  interval  (u  ,v  ) ,  R  is  a  maximal  effort  curve.  Because 

n  n 

♦  ♦  —  4  4-  4-  -v 

P  >  y(P)  -  EP  >  0  whenever  U  <  P  <  V  ,  P  must  be  strictly  increasing  on  [u  ,v  ] 

n  n  n  n 

until  the  time  at  which  P  reaches  VR.  Of  course,  if  n  *  N,  P  may  never  reach  VN,  in 
which  case  P  is  strictly  increasing  on  [u^,®).  Let  C^,  Un  -  <  Cj  < •  •  •  <  CK  *  Vn , 

denote  the  extrema  of  y  on  [Un,Vn) ,  so  that  K  is  odd,  and 

for  k  even,  y  is  strictly  increasing  on  , 


(42) 


strictly  decreasing  on  [0^,0^^] 


Let  t,  be  the  time  at  which  R(t,  )  *  C.  .  Further,  let  Q  be  the  function  obtained  by 


stretching  P  at  each  for  odd 

k  <  K,  and  then  restricting  the  result 

to  [u  ,v  ] .  The  amounts  of  the  stretches 
n  n 

are  uniquely  determined  by  the  requirement 
Q(t^)  =  R(t^)(=  C^)  for  k  even.  It 
is  clear  from  this  construction  and  (27) 
that,  for  k  even, 


r 

k-1 


<  Q  t  R  < 


on 


[t, 


-,'V 


<  R  <  Q  C 


k+  1 


-Vi’ 


Thp»se  i neaiia  1  i t  i *><;  and  r42)  im*' 1  y  that 

y ( R ( t ) )  >  v  ( Q  ( t )  )  on  (u  ,v  1 

n  n 


(43) 


The  integral  of  Y(P)  over  the  times  at  which  P  is  strictly  increasing  is  equal  to  the 


analogous  integral  for  y(Q).  Thus,  since  y(Ct)  >  y(C  ), 

*  K 

V 

Jj*  lY<Q(t>)  -  y  (P(  t) ))  dt  >  0  , 

u 

n 

and  we  conclude  from  (43)  that 

v 

Jj1  [y(R( t) )  -  y (P(  t) )  ]  dt  >  0  , 

u 

n 

for  n  <  N.  Thus  to  show  that  11m  Inf  I(J,p)  >0  we  have  only  to  show  that 

lim  inf  (y(Q(t))  -  y(£(t))]dt  >  0 

T  *  u 

n 

To  verify  (45)  let  (k  odd)  be  the  length  of  the  stretch  in  Q  at  and  let  L 

denote  the  sum  of  the  i Then  the  integral  in  (45)  is  equal  to 

t-l  [s  "  y(p<t))]<3t  +  t  \ty(ck)  -  s]  . 

odd  k 

•* 

The  integral  tend3  to  zero  as  T  -►  «,  since  P(t)  ♦  Pg  as  t  ♦  ”,  while  the  sum  is 
>  0,  since  y(C^)  *  S.  Thus  (45)  holds. 

The  remainder  of  the  proof  of  (41),  that  lim  inf  I(JT)  >  0  when  R  f  P  on  J, 
safely  be  omitted.  This  completes  our  proof  of  Assertion  4. 

We  now  return  to  the  proof  of  (38)  and  consider  separately  two  cases. 

*  > 

Case  1  (R  *  P  ).  Here,  by  (37),  R  +  P,  and  (38)  follows  from  Assertion  4  with 
R  =  P  . 

*  -* 

Case  2  (R  y  P  ) .  Choose  T  >  ( >  u^)  and  consider  the  integral 

I0  »  /g  (y{ P* ( t) )  -  v(R(t)))dt  . 


(44) 


(45) 


can 
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The  inteqral  of  y(R)  over  the  free-qrowth  portions  of  R  is  equal  to  the  analoqous 


* 

inteqral  for  P  •  Similarly,  the  inteqral  of  y(R)  over  the  maximal-effort  seqments  of 
R  on  [0,11^]  equals  its  counterpart  for  P  on  {0,uNl.  Thus  7q  =  ^2  “  1 3 •  where 

J1  =  /I  y(R(t))rtt,  12  =  fl  y  (P*(  t ) )  dt 
UN  UN 

and  Is  the  inteqral  of  y(R)  over  those  subintervals  of  I  ( cf .  (39)^)  on  which  P 

is  constant.  The  measure  of  these  sub-intervals  is 

T  =  %  •  UN  • 

Thus,  since  y(R(t))  <  S  on  I,  we  have  the  bound 

7  3  <  tS  -  A 

with 


(1) 

A  =  0 

when 

y(R(t) ) 

=  s 

on 

all  of  the  above  subintervals 

(2) 

A  >  0 

when 

y  ( R  ( 1 1  ) 

<  s 

on 

at  least  one  such  subinterval 

Finally, 

since 

R(t+x) 

=  P*(t) 

for 

t 

>  UN ' 

I2  -  =  £_T  y(P#(t)  )dt  . 


In  view  of  the  above  remarks. 


Y  ( p  ) 
T 


Y  (R)  >  P  (T-t)  -  P  (T) 
T 


|  !y (P  (t) )  -  S) dt  t  A 


* 

and,  since  P  (t)  ♦  P  . 

S 

lim  inf  (Y  (P*)  -  Y  (R)l  >  A  ?  0 
T  T 

T  ♦  ® 


UFO 


For  A  >  0  this  result  and  Assertion  4  imply  (38).  Thus  assume  A  *  0.  Since  R  ¥  P# 
there  is  at  least  one  interval  ft  in  I  with  R  constant  at  S,  so  y(R)  -  ER  * 

S  -  ER  >  0  on  ft.  Since  P  *  P  e  A  on  I  and  R  fails  to  satisfy  (A3)  on  ft  c  I, 
p  /  R  on  ft.  Thus  (38)  follows  from  (46)  and  Assertion  4.  This  completes  the  proof, 
d.  Discounting. 

The  usual  method  of  treating  an  infinite  horizon  problem  is  to  introduce  a  discount 
factor  e  ***  into  the  yield.  The  total  yield 

J*  e_<tE(t)P(t)dt 

is  then  defined  for  all  bounded  paths  P  and  may  be  written  as 

YIP)  =  /”  e_,Ctw(P(t))dt  +  PQ 


where 


w (P)  =  y(P)  -  <P  .  (47) 

The  corresponding  optimal  harvesting  problem  for  unconstrained  effort  consists  in 
maximizing  Y  over  the  class  of  all  bounded  paths  which  are  consistent  with  (Al)  and 
(A2).  If  w  is  a  C1  function  on  [0,®)  with  a  strict  global  maximum  at  PM  >  0  and  if 

w  has  no  other  extrema  on  (0/®)  (48) 

then  the  solution  of  this  problem  is  analogous  to  the  solution  of  the  nondiscounted  problem 
of  section  b.  ( cf . ,  e.g.,  Clark  and  Munro  [1]  and  Clark  [2],  chapter  2.)  We  now  show,  by 
example,  that  if  the  rather  restrictive  assumption  (48)  is  eliminated,  the  optimal  path 
not  be  a  most  rapid  approach  to  P^. 
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P  P  e 

0  0 


Figure  15.  The  function  w. 


Let  Pg  >  0,  x  >  0,  and  M  >  0  be  fixed.  Further,  let  e  >  0  be  arbitrarily  small, 
let  PM  =  PM(e)  =  p0/e,  and  let  w  “  we  be  a  c1  function  on  [0,°»)  with  the  following 
properties:  w  has  a  strict  global  maximum  at  PM,  w  >  0  on  [Pg,pMl, 

w(PM)  »  M,  w(PQ)  -  M  -  e  ,  (49) 


M  -  e  on  [PQ,PQe  1 


'  2  2 

€  -€ 

V  e  on  (P  e  ,  P  e  ] 

U  M 


we  will  show  that  for  e  sufficiently  small  the  path  P  =  is  better  than  the  path  P 

defined  in  ( 7 ) . 

To  see  this  let  Z  =  Z_  he  the  free-growth  curve  starting  from  Pn  at  t  =  0,  so 


that ,  by  ( 47 ) , 


Z  =  w(Z)  +  <Z  . 


Let  tM  =  tM(e)  be  the  time  at  which  Z  reaches  PM  ( cf .  (7)).  Further,  let  tg  *  tQ(e) 

P2  _p2 

and  t^  =  t^(0)  denote  the  times  at  which  Z  reaches  P^e  and  P^e  ,  respectively, 

so  that  tg  <  tj  <  tM.  Assume  that  @  <  tc.  Then  (SI)  yields 


PZ  <  7.  <  M  >  *7,  on  [n,t  1 


-e2 

Since  M  and  <  are  Independent  of  @,  while  P^e  becomes  infinite  as  S  ♦  0,  we 
may  conclude  from  the  second  inequality  in  (52)  that 

t(  *  *  is  e  <  d  ,  (">3) 

et0  e2 

Similarly,  the  first  inequality  in  (52)  yields  Z(tg)  >  P^e  ,  and,  since  Z( t^ )  «  PQe  , 

it  r allows  that  t  5  e;  thus 
0 


tQ  ♦  0  as  e  ♦  0  .  (54) 

In  view  of  (50),  w(Z(t>)  is  bounded  above  by  M  -  e  for  0  <  t  <  t  ,  by  e  for 

t.  <  t  <  t  ,  and  by  M  for  t  <  t  <  t  t  thus,  using  (7)  and  (49)., 

0  1  1  M  ^ 


**t« 


Y(P*)  -  Y  (P )  <  /  1  (2e-M)e_<tdt  +  C  ee_,Ctdt 

fc0 


( 26-M) e 


+  (M-3e)e 


By  (S3)  and  (54)  the  right  side  of  the  above  relation  has  the  limit  -M/<  as  6  ♦  0, 

n  * 

Hence  we  have  only  to  choose  e  sufficiently  small  to  insure  Y(P)  >  Y(P  ). 


3.  The  age-dependent  theory  without  harvesting, 
a.  Basic  equations.  Persistent  age  distributions. 

We  work  within  the  framework  of  the  nonlinear  theory  developed  in  [5).  This  theory  is 
based  on  the  equations 


+  M^ti+  u(a,P(t))p(a,t)  ,  o  , 

da  dt 

P(  t )  =  f°°  p(  a ,  t)  da  , 


B(t)  =  p(0,t)  *  /  B(a,P(t) ) p(a,t)da  , 


(55) 


where 


p(a,t)  is  the  age  distribution  (the  number  of  individuals,  per  unit  age, 
of  age  a  >  0  at  time  t  >  0)j 
P(t)  is  the  total  population; 

B(t)  is  the  birth-rate; 


-2R- 


V  is  the  death  function  (y(a,P)da  is  the  probability  of  dying  in 


the  age  interval  (a,a+da)  when  the  population  is  P); 

0  is  the  birth  function  (0(a,P)  is  the  expected  number  of  births 

to  a  person  of  age  a,  per  unit  time,  when  the  population  is  P) . 

Here,  however,  we  assume  that  dependence  on  total  population  is  confined  to  the  death 
function*  More  specifically,  we  assume  that1 


B(a,P)  *  0(a)  , 

u(a,P)  »  u  (a)  +  u  (P) 
n  e 


(56) 


In  (56).,,  y  (a)da  represents  the  probability  of  dyinq  of  natural  causes  durinq  (a,a+da), 
z  n  - 

while  y  (P)da  is  the  probability  of  death  due  to  environmental  factors  (crowdinq,  etc.) 
e  - 

during  the  same  interval.  In  view  of  this  interpretation. 


n(a)  =  exp{^  Un  ( a)  da} 

is  the  probability  of  living  to  age  a  in  a  perfect  environment?  that  is,  an  environment 
with 

Vi  =  0  .  (57) 

e 

Assumption  (56 ^  is  special,  as  the  effects  of  the  environment  are  independent  of  age?  it 

might  be  appropriate,  for  example,  to  a  population  in  the  presence  of  predators  which 

indiscr iminatly  eat  individuals  of  all  aqes. 

We  assume  y  ,  y  ,  and  0  are  C1  functions  on  [0,«>)  with  0  e  l\o,«)  and 
n  e 

J*  *(a) 0( a) da  >  1  .  (5B) 


1 

These  assumptions  with  u  =  0  and  0(a)  a  sum  of  terms  of  the  form  ba^e”  ia  were 
utilized  by  Ourtin  and  MacCamy  [5,H]  to  reduce  the  system  (55)  to  ordinary  differential 
equations.  The  more  general  form  (56)  was  introduced  by  Coleman  and  Simmes  [41- 
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The  left  side  of  (58)  represents  the  net  reproduction  rat-e1 2 3  in  a  perfect  environment.  A 
consequence  of  (58)  is  that  in  such  an  environment  the  population  ultimately  grows.  This 
assumption  therefore  yields  a  population  which  is  amenable  to  harvesting. 

Because  of  (58),  the  equation 


/“  Ti(a)6(a)e-rada  =  1 


has  exactly  one  real  solution  r,  and  r  >  0.  We  call  r  the  natural  growth  rate.  For  a 
perfect  environment  the  function 

p(a,t)  =  CTi(a)e  r*P(t)  (59) 


with 


P(t)  =  P(0 )ert 


and  C  chosen  so  that 


C  fg  ir(a)e  ra<3a  “  1  (60) 

is  a  solution  of  (55)  ((56)  and  (57)).  We  call  such  solutions  persistent  age 
distributions.  Their  importance  is  that  they  represent  the  asymptotic  behavior  of  general 

2 

solutions  for  large  time:  given  any  solution  p  there  is  a  constant  Cg  such  that,  for 
each  a , 

-ra  rt 

p(a,t)  -  Cqtt(  a^e  6 
3 

as  t  ♦  «. 


1 

Cf.,  e.g.,  Key  fit?.  [9],  p.  102. 

2 

Fo 1 1 er  M  n }  . 

3 

f(M  -  «( t )  is  t  -►  «  signifies  that  f(t)/g(t)  +  1  as  t  -*•  ». 
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Once  8  is  known  the  total  population 


P(t>  -  /"  ?(a,t)d a 


(65) 


corresponding  to  5  is  easily  calculated  using  the  relations 

P(t)  =  /Jj  ir(a ) 8 ! t-a ) da  +  Yft)  , 

¥(t)  =  r< a-t,a)^ (a-t)  da  . 


(66) 


It  is  important  to  emphasize  that  £,  8,  and  P  are  the  age  distribution,  birth-rate,  and 
total  population  that  would  prevail  in  a  perfect  environment. 

To  derive  (64)  and  (66)  from  (63)  we  integrate  (63)1  along  characteristics  (t  =  a  + 
constant).  The  result. 


£(a,t)  » 


{it(  a)8(  t-a) ,  0  <  a  <  t 

ir(  a-t,  (  a-t)  t  4  a 


(67) 


when  combined  with  ( 63 ) 2  and  (65),  yields  (64)  and  (66). 

Next,  we  integrate  (62)  with  respect  to  age  and  use  < 5 5 ) 2  and  (65)  to  conclude  that 


P(t)  =  P(t)exp{/*  ue(PU))dX) 


(68) 


If  we  differentiate  this  relation  with  respect  to  t  and  divide  the  result  by  P(t),  we 

find  that  P  satisfies  the  ordinary  differential  equation 

P  =  (y  -  U  (P)]P  .  (69) 

e 


Pit) 

Pit) 


Here 

y(  t) 

represents  the  instantaneous  growth  rate  in  a  perfect  environment.  Of  course,  P 
satisfies  the  initial  condition 


(70  ) 
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L  v  (a) da 


(71) 


P(0) 


The  above  analysis  yields  the  following  procedure  for  solving  the  initial-value 
problem  (55),  (56),  and  (61): 

(a)  Solve  the  linear  inteqral  equation  (64)  for  8(t). 

(b)  Compute  P(t)  from  (66)  and  >(t)  from  (70). 

(c)  Solve  the  differential  equation  (69)  -  with  initial  condition  (71)  -  for  P(t). 

(d)  Compute  !j(a,t)  from  (67)  and  p(a,t)  from  (62). 

4.  Harvesting  of  age-structured  populations, 
a.  Basic  equations.  The  optimal  harvesting  problem. 

We  consider  harvesting  with  effort  Eft)  independent  of  aqe.1  Thus  individuals  of 
age  a  are  harvested  at  a  rate  E(t)p(a,t),  per  unit  age  and  time,  and,  in  place  of 
( 55  )  g ,  we  have  the  balance  equation 

Miill  +  +  ( p(a,P(t) )  +  E(t))p(a,t)  =  0  .  (72) 

The  underlying  equations  are  therefore  (72)  and  (55)2  ^  supplemented  by  the  initial 
condition  (6t),  (We  continue  to  assume  that  u  and  ft  satisfy  (56).) 

We  will  use  the  procedure  discussed  in  Section  3  to  reformulate  the  problem.  We 
therefore  define 


FJa  ,t)  =  p(a,t)exp{/^  (u  (P(Jt))  +  EU)]dX)  .  (71) 

Then  £  satisfies  the  linear  system  (63),  6(t)  =  F,(0,t)  is  again  the  solution  of  the 
linear  integral  equation  (64),  and  P(t),  defined  by  (65),  is  again  given  by  (66).  It  is 


1 

Cf.  Remark  3. 
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important  to  note  that  £,  8,  and  P  represent  the  age  distribution,  birth-rate,  and 


total  population  that  would  prevail  in  the  absence  of  harvesting  and  environmental  effects. 

As  before,  the  behavior  of  the  actual  population  P(t)  is  governed  by  an  ordinary 
differential  equation.  Indeed,  the  exact  same  steps  used  to  derive  (69)  now  lead  to 

P  =  y(t,P>  -  EP  (74) 

with 

y(t,P)  =  (y(t)  -  p  (P)]P  .  (75) 

e 

In  view  of  (74),  the  yield 

£  /“  E(t) p(a,t) dadt  -  G  E(t)P(t)dt  (76) 


on  a  time  interval  (0,T)  is  given  by  the  functional 

VP)  =  p0  ■  p<T_>  +  G  y(t«p(t>>dt  •  (77> 

Roughly  speakinq,  our  optimal  harvesting  problem  consists  in  maximizing  (77)  sub jet  to 
certain  constraints.  This  problem  can  be  attacked  using  the  following  procedure: 

(a)  Solve  the  linear  integral  equation  (64)  for  B(t). 

(b)  Compute  P(t)  from  (66)  and  y(t)  from  (70). 

(c)  Solve  the  optimal  control  problem  defined  by  the  differential  equation  (74) 
and  the  functional  (76). 

Unfortunately,  because  of  the  dependence  of  y(t,P)  on  t,  the  optimization  problem 
of  (c)  doe^  not  fit  within  the  framework  discussed  in  Section  1.  There  is,  however,  an 
important  class  of  problems  for  which  this  dependence  disappears.  This  class  corresponds 
to  initially  persistent  aqe  distributions ,  and  we  shall  study  it  in  the  next  section, 
b.  optimal  harvesting  for  initially  persistent  age  distributions. 

As  noted  at  the  end  of  Section  3a ,  when  the  environment  is  perfect  and  harvesting 
absent,  fV*  vip  distribution  r  is  ultimately  persistent.  Therefore,  if  the  population  has 
been  evnlviner  over  a  Iona  period  of  time,  and  if  all  past  harvesting  has  been  with  effort 
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independent  of  age,  then  it  seems  reasonable  to  assume1  that  initially  p  will  have  the 
age  structure  indicated  in  (59).  We  therefore  assume  that  the  population  is  Initially 
persistent  in  the  sense  that 

(a)  »  CPgit(  a)e  ra  , 

where  C  and  Pg  are  constants  with  C  chosen  so  that  (60)  holds,  while  r  is  the 
natural  growth  rate.  The  initial-value  problem  (63)  then  has  the  unique  solution 


£<a,t)  -  Cir(a)e-ra  R  t) 
P(t)  -  PQert  , 


(78) 


so  that,  by  (70), 


y( t)  -  r  .  (79) 

Further,  the  counterpart  of  (68)  in  the  present  circumstances  is  simply  (68)  with  E( X) 

added  to  iMP(X));  thus  (78)1  and  (73)  imply  that 

p( a , t)  -  Cit(a)e  rap(t)  .  (80) 

In  view  of  (75)  and  (79),  y(t,P)  is  independent  of  ts 

y(P)  -  (r  -  p  (P))P  .  (81) 

e 

Thus  the  differential  equation  (74)  reduces  to 

P  -  y(P)  -  EP  ,  (82) 


while  the  functional  (77)  takes  the  form 


VP)  =  PQ  -  P(T~ )  +  y(P(t))dt 


(83) 


Comparing  (82)  and  (83)  with  (3)  and  (4),  we  see  that  for  an  initially  persistent 
population  the  optimal  harvesting  problem  reduces  to  the  (age-independent)  problem 
discussed  in  Section  2.  Therefore,  if  y  satisfies  the  hypotheses  of  Section  2,  the 
theorems  of  Sections  2b  and  2c  hold. 


1 

Cf.  (73). 
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Remark  3.  Our  results  extend  trivially  to  situations  in  which  individuals  of  age  a 


have  an  economic  value  g(a),  where  g  is  a  nonnegative,  L°  function  on  (0,“).  In 
this  case  the  yield  (76)  is  replaced  by 


f0  /”  g(a)E(t)p(a,t)dadt 


f1  E(t)G(t)dt 


(84) 


where 


G(t)  -  g(a) p(a,t)da  , 


and  (80)  implies  that 


G(t)  =  (C  f“  g(  a)  ir(  a)  e  r*da)P(t)  =■  C,P(t) 

with  C,  >  0  constant.  Thus  (84)  is  a  constant  times  the  yield  (76),  and  the  corresponding 
optimisation  problem  reduces  to  that  already  discussed.  When  the  age  distribution  is  not 
initially  persistent  this  reduction  does  not  take  place;  the  resulting  problem  is  treated 
in  (11]. 
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